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FOREWORD 

A combined  numerical-experimental  study  was  performed  in  the  Turbo- 
Structures  Research  Lab  (TSRL)  between  1 June  1977  and  1 December  1977. 

The  purpose  was  to  Investigate  the  vibratory  response  of  an  Internally 
clamped  annular  disk  under  both  static  and  rotational  conditions. 

A 10  Inch  diameter  steel  disk  with  an  Integral  hub  was  manufactured. 
Holographic  Interferometry  was  used  to  experimentally  determine  the  fre- 
quencies and  mode  shapes  for  the  first  five  diametrical  modes  and  one 
complex  mode  under  nonrotating  conditions.  Incorporating  an  image  derota- 
tor and  a double-pulsed  ruby  laser,  resonant  conditions  and  mode  shapes 
were  again  determined  for  the  above  modes,  but  this  time  the  interferograms 
were  constructed  while  the  disk  was  rotating  at  specified  speeds. 


The  finite  element  computer  program  NASTRAN  was  used  to  model  the 
steel  disk  and  analytically  determine  the  mode  shapes  and  frequencies. 

For  the  nonrotating  case,  all  mode  shapes  and  frequencies  were  determined 
and  compared  to  the  corresponding  modes  found  experimentally.  Under  condi- 
tions of  rotation,  only  the  third  diametrical  mode  was  studied  as  a function 
of  centrifugal  loading  because  of  the  large  requirement  for  computer  central 
memory  and  associated  expense.  The  resulting  frequency  predictions  were 
compared  to  the  experimental  findings  with  a deviation  of  about  4%. 

The  work  was  performed  In  the  Turbine  Engine  Division  of  the  Air  Force 
Aero- Propul slon  Laboratory,  Air  Force  Systems  Command,  Wrlght-Patterson  Air 
Force  Base,  Ohio,  under  Project  3066,  Task  12,  and  Work  Unit  21.  The  effort 
was  conducted  by  John  E.  Horner,  Dr.  James  C.  MacBaln, and  William  A.  Stange 
of  the  Propulsion  Branch. 

The  authors  are  Indebted  to  Mr.  Bruce  Tavner  for  his  competent  technical 
assistance  In  the  laboratory  and  to  Miss  Helen  Davis  for  typing  the  manuscript. 


Ill 


II  DISCUSSION  4 


II  DISCUSSION  4 

1.  Disk  Vibration  4 

2.  Experimental  Procedures  - Static  Conditions  6 

3.  Experimental  Procedures  - Rotating  Conditions  9 

4.  Alignment  13 

5.  Numerical  Analysis  14 

III  RESULTS  17 


1.  Experimental  Results 

17 

2.  Finite 

Element  Analysis  and  Theoretical  Correlation 

21 

CONCLUSIONS 

24 

APPENDIX  A 

Burst  Speed  Calculation 

27 

APPENDIX  B 

Hub  Effects  on  Frequencies 

28 

APPENDIX  C 

Equipment 

30 

REFERENCES 

32 

v 


BMMlfci 


AFAPL-TR-78-62 


LIST  OF  ILLUSTRATIONS 


FIGURE  PAGE 

1 Frequency  Speed  Diagram  35 

2 Experimental  Set-Up  36 

3 Schematic  of  Set-Up  36 

4 Experimental  Set-Up  for  Image  Derotated  Holographic 

Interferometry  3/ 

5 Image  Derotator  38 

6 Lathe  Motor  and  Disk  38 

7 Folded  Abbe  Inverting  Prism  39 

8 Steel  Disk  40 

9 Vibration  Instrumentation  40 

10  Resonant  Mode  Shapes  In,  2n,  3n  for  the  Stationary 

Disk  41 

11  Resonant  Mode  Shapes  4n,  5n,  and  (2n,  Is)  for  Sta- 
tionary Disk  42 

12  Resonant  Excitation  Frequency  as  a Function  of  Disk 

Speed,  Modes  In,  2n,  and  3n  43 

13  Resonant  Excitation  Frequency  as  a Function  of  Disk 

Speed,  Mode  4n  and  5n  44 

14  Double  Exposure  Holograms  of  Disk  Modes  In,  2n,  and 

3n  at  Specified  Excitation  Frequency  and  Disk  Speed  45 


vl 


AFAPL-TR-78-62 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


FIGURE 


15  Double  Exposure  Holograms  of  Disk  Modes  4n,  5n,  and 
(2n,  Is)  at  Specified  Excitation  Frequency  and 
Disk  Speed 

16  Double  Exposure  Hologram  of  4n  Disk  Mode  at  High 
RPM 

17  Finite  Element  Model 


18  Finite  Element  Resonant  Mode  Shapes  In,  2n,  and 
3n 

19  Finite  Element  Resonant  Mode  Shapes  4n,  5n,  and 
(2n,  Is) 

20  Disk  Natural  Frequency  as  a Function  of  Speed 

21  Complex  Mode  Shape  of  Stationary  Bladed  Disk 


46 

47 

48 

49 

50 

51 

52 


22  Clamping  Fixture 


53 


vl  1 


AFAPL-TR-78-62 

LIST  OF  TABLES 

TABLE  PAGE 

1 Experimental/Analytical  Correlation  - Steel  Disk  33 

] 

2 Experimental/Analytical  Correlation  - Aluminum  Disk  34 


AFAPL-TR-78-62 


\ 


n 

s 

U 

P 

Fr 


t 

A 

6. 


n 

e 

a 

b 

T 


U) 


P 

E 

v 

D 

6 


NOMENCLATURE 


Nodal  Diameter 
Nodal  Circle 

Out-of-Plane  Displacement 
Excitation  Force 
Excitation  Frequency 

Time 

Resonant  Amplitude 

Absolute  Angular  Coordinate  of  Disk 

Disk  Vibration  Frequency 

Angular  Velocity  of  Disk 
Angular  Coordinate  of  Disk 
Disk  Outer  Radius 


Disk  Inner  Radius 
Thickness 


Eigenvalue 
Frequency 
Mass/ Area 
Young's  Modulus 
Poisson's  Ratio 
Flexural  Rigidity 
Density 

Average  Stress 
Critical  Speed 


lx 


I 


4 


AFAPL-TR-78-62 


SECTION  I 
INTRODUCTION 

Rotating  disk-like  structures  are  encountered  throughout  the 
engineering  development  of  advanced  gas  turbine  engines  and  have  been 
the  object  of  study  since  the  appearance  of  the  steam  turbine  (Reference 
1).  They  possess  several  families  of  modes  of  free  vibration  comprised 
of  both  circumferential  and  diametral  modes,  each  having  an  Infinite 
number  of  possible  resonant  frequencies.  Each  resonant  response  corres- 
ponds  to  a deflection  pattern  comprised  of  nodal  lines  which  can  rotate 
relative  to  the  disk  or  remain  fixed  to  the  disk  depending  on  the  rela- 
tive phase  shifts  as  the  structure  passes  through  resonance.  These 
responses  are  difficult  to  characterize  because  the  maximum  stress  or 
strain  location  is  constantly  changing  for  each  mode  of  vibration  In  a 
particular  family. 

Resonant  frequencies  of  the  disk  are  affected  by  the  various 
operating  conditions  such  as  rotational  speed,  temperature,  pressure, 
and  boundary  effects.  Fatigue  damage  can  accumulate  quite  rapidly  when 
the  structure  operates  within  a resonant  band,  especially  when  It  does 
not  possess  sufficient  damping  or  when  It  passes  too  slowly  through  this 
condition.  The  excitation  forces  which  can  excite  resonance  in  a rotat- 
ing disk  are  numerous.  Dynamic  unbalance,  oil  pumps,  gear  meshing,  ball 
bearings,  engine  struts,  pressure  differentials,  and  a variety  of  other 
sources  may  excite  a disk  at  frequencies  from  fractional  to  multiple 
orders  of  engine  speed.  They  can  excite  disk-like  structures  Into  large 
deflections  where  nonlinearities  can  occur.  Under  dynamic  conditions, 
the  nonlinear  systems  may  respond  to  frequencies  equal  to  and/or  In 
multiples  of  the  forcing  frequency.  Their  occurrences  are  accompanied 
by  fatigue  damage  accumulating  at  a rate  higher  than  expected,  and  at  a 
resonance  which  may  occur  at  frequencies  removed  from  the  forcing  fre- 
quencies. 
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Analytically,  it  Is  difficult  If  not  Impossible  to  currently  pre- 
dict the  structural  response  of  disk-like  structures  operating  vithln 
a turbine  engine  environment.  In  order  to  establish  and  verify  the 
vibrational  characteristics  of  these  structures,  they  must  be  manufactured 
and  tested  experimentally.  This  is  usually  accomplished  by  instrumenting 
with  strain  gages  and  recording  the  output  via  telemetry  or  slip  rings. 

This  method  suffers  from  mechanical/electrical  linkage  problems  and,  per- 
haps more  often,  from  inadequate  placement  of  the  gages  themselves.  Data 
reduction  and  interpretation  of  results  are  difficult  and  cumbersome. 

Through  the  advent  of  holography,  alternative  methods  have  been 
established  to  determine  the  response  of  structures  as  a function  of 
various  forcing  frequencies  under  nonrotating  conditions.  Mode  shapes 
and  frequencies  are  readily  determined  for  the  full  surface  of  the  struc- 
ture, a feature  which  often  leads  to  quick  identification  of  areas  of  high 
stress  concentrations.  However,  effects  of  changing  boundary  conditions  on 
structural  stiffness,  and  the  resulting  resonant  conditions  as  a function 
of  centrifugal  loading  could  not  be  investigated  because  of  the  rigid  re- 
quirements of  no  relative  movement  required  to  conduct  holographic  inter- 
ferometry. 

As  early  as  1970  (Reference  2),  researchers  have  attempted  to  apply 
holographic  interferometry  to  rotating  objects.  Success  was  met  through 
an  Air  Force  contract  granted  to  United  Technologies  Research  Center  in 
developing  an  optical  system  for  dynamic  analysis  of  rotating  structures 
(Reference  3).  The  method  uses  an  image  derotator  to  optically  compensate 
or  stop  the  motion  of  the  object,  giving  it  the  appearance  of  a stationary 
object.  It  gives  continuous  motion  compensation  of  the  object  and  allows 
double  pulse  holography  to  be  performed  while  the  object  is  rotating  and 
viewed  through  the  image  derotator.  Once  the  holographic  film  is  developed, 
the  out-of-plane  displacement  and  resulting  mode  shapes  can  be  photographed 
and  studied  when  the  film  is  illuminated  by  a laser  light.  The  system  was 
delivered  to  the  Air  Force  and  installed  in  the  Turbo-Structures  Research 
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Lab  (TSRL).  The  present  work  was  to  Install  the  system,  experimentally 
determine  the  resonant  mode  shapes  and  frequencies  for  a rotating  disk, 
and  verify  these  results  through  the  utilization  of  numerical  techniques. 
Discussed  will  be  the  theoretical  aspects  of  disk  vibration,  mode  shape 
determinations  under  nonrotating  conditions,  the  same  modes  as  a function 
of  Increasing  speed  or  centrifugal  loading  and  the  results  of  the  finite 
element  program  NASTRAN  used  to  verify  the  experimental  results. 
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SECTION  II 
DISCUSSION 


1.  DISK  VIBRATION 

Assume  a circular  disk  fixed  at  Its  center  and  unconstrained  along 
Its  outer  periphery.  The  modes  of  vibration  may  be  composed  of  nodal 
diameters  (n),  nodal  circles  (s),  or  combinations  of  nodal  diameters  and 
circles  (Reference  4).  Discussion  shall  concern  those  modes  containing 
only  nodal  diameters,  though  arguments  are  equally  applicable  to  a com- 
bination of  nodal  diameters  and  circles.  Modes  comprised  of  nodal  circles 
shall  not  be  addressed  In  this  discussion.  They  require  axlsymmetrlc 
forcing  functions  to  drive  them  which  are  rarely  encountered  In  practice. 

The  vibration  of  the  nonrotating  disk  vibrating  In  the  nth  nodal 
diameter  excited  by  a stationary  force  varying  sinusoidally  In  time, 

Pcos  FEt,  can  be  described  by: 

W = A sin  n 0Q  cos  pnt  (1) 

where  W Is  the  transverse  displacement,  and  A sin  n 0Q  indicates  the 
shape  of  vibration;  A Is  the  resonant  amplitude,  which  Is  a function 
of  the  disk  radius  and  proportional  to  the  magnitude  of  the  excitation 
force,  P.  The  absolute  angular  coordinate  of  the  disk  relative  to  Iner- 
tial space  Is  denoted  by  0Q.  (See  Figure  1 for  clarification  of  notation.) 
Likewise,  cos  pnt  Is  the  harmonic  movement,  with  pn  being  the  vibration 

frequency  and  t the  time.  Resonance  will  occur  when  the  excitation  fre- 
quency, Fe,  Is  equal  to  pn> 

We  can  use  a trigonometric  Identity  to  put  Equation  1 In  the  form: 

W = A sin  (n  eQ  - pnt)  + A sin  (n  0Q  + pnt)  (2j 
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which  can  be  considered  as  the  superposition  of  two  traveling  waves 
of  shape  A sin  n0Q.  Now,  If  the  disk  Is  rotated  at  an  angular  velocity, 

n,  the  absolute  angular  coordinate,  60,  Is  related  to  the  angular  position 
of  the  disk,  0,  by  0Q  = 0-flt.  Substituting,  Equation  2 becomes: 

W = A sin  [n0  - (pp  + nfl)  t]  + A sin  [n0  + (pn  - nfi)  t]  (3) 

Looking  at  Equation  3,  we  see  that,  under  conditions  of  rotation,  the  disk 
resonance  will  occur  when  the  excitation  frequency,  FE«  Is  such  that  F£  = 

Pn  + n Q.  This  Is  saying  that  two  excitation  frequencies  will  excite 

the  same  resonant  mode  at  a given  rpm.  For  example,  the  2n  mode  of  vibra- 
tion will  be  excited  Into  resonance  by  two  different  excitation  frequencies, 
one  that  Is  2 ft  greater  than  the  2n  natural  frequency  and  one  that  Is  2ft 
less  than  the  2n  natural  frequency.  A typical  frequency  speed  diagram  demon 
stratlng  this  phenomena  Is  shown  In  Figure  1. 

Here,  line  AB  represents  the  natural  frequency  of  the  disk  for  the 
nth  mode  of  vibration  as  a function  of  speed,  n.  The  Increase  In  fre- 
quency Is  attributed  to  the  centrifugal  stiffening  as  a result  of  rotation. 
Line*  AC  and  AD  represent  the  upper  and  lower  resonant  conditions, 

Fe2  = Pn  + nfiQ  and  F£1  = Pn  - nflQ  respectively.  Thus,  at  speed 

resonance  will  exist  as  Indicated  by  the  applied  frequencies  FE2  and  F^. 

Experimentally,  line  AB  Is  easily  established  by  bisecting  the  upper  and 
lower  curves  of  AC  and  AD  which  are  determined  by  monitoring  disk  trans- 
verse displacement  for  Increasing  values  of  rpm.  This  line,  which  depicts 
the  natural  frequency  of  the  disk,  was  later  verified  for  the  3n  mode  by 
utilizing  finite  element  techniques.  At  point  D,  when  n«c  * pn>  an 

Interesting  phenomenon  occurs.  Here  the  speed  of  rotation  equals  fi  and 


5 


AFAPL-TR-78-62 


resonance  can  result  from  a stationary  force  P,  which  has  an  excitation 
frequency,  FE,  equal  to  zero,  (This  Is  equivalent  to  a stationary  force 
P traveling  at  nil  around  the  rim  of  a stationary  disk.)  In  a turbine 
engine,  this  condition  can  occur  quite  readily  within  the  speed  range  of 
a given  compressor  stage.  Attempts  to  operate  In  this  region  for  pro- 
longed periods  of  time  are  therefore  avoided.  For  speeds  In  excess  of 
fic,  two  more  resonant  frequencies  are  found  for  each  value  of  speed. 

Line  DE  Is  the  mirror  Image  of  the  extension  of  AD  and  represents  the 
case  when  Pn  = - (F£  - nft),  which  provides  a positive  value  for  pn  when 
Fe  < nfi. 

2.  EXPERIMENTAL  PROCEDURES  - STATIC  CONDITIONS 

A circular  disk  was  manufactured  from  1018  cold-rolled  steel.  It 
had  a constant  thickness  of  0.0625  Inches,  with  an  outside  diameter  of 
10  Inches.  Integrally  formed  with  the  disk  was  a center  hub  which  was 
1.0  Inch  In  diameter  and  1.0  Inch  thick,  and  had  a center  bore  hole  of 
0.5  Inches  for  sliding  onto  the  shaft  of  a drive  spindle.  The  transition 
from  the  hub  to  the  disk  was  made  with  a 1/16  Inch  radius  fillet. 

Initially,  the  disk  was  statically  mounted  to  an  Isolation  table 
and  Its  first  five  diametral  modes  of  vibration  (In,  0s)  through  (5n,  0s) 
and  the  second  coupled  mode  of  vibration  (2n,  Is)  determined  using  time 
average  holography.  In  classifying  the  different  modes  of  vibration,  n 
refers  to  the  number  of  nodal  diameters  and  s refers  to  the  number  of 
nodal  circles,  not  Including  the  Inner  fixed  boundary. 


Under  static  conditions,  holography  Is  an  Ideal  method  of  experi- 
mentally determining  the  natural  modes  and  frequencies.  As  early  as 
1965  (Reference  5),  holographic  Interferometry  was  first  used  to  study 
vibrating  objects.  It  has  since  proved  useful  throughout  the  aerospace 
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Industry  for  vibration  analysis  of  critical  components  such  as  compressor 
blades  and  disks.  Unlike  strain  gages  and  photoelastic  coatings,  holo- 
graphy Is  a non-lntruslve  optical  technique  that  does  not  affect  the  mode 
shapes  or  response  of  the  structure.  The  object  can  be  of  any  shape  with 
the  final  results  being  a full  field  picture  of  the  vibrating  displace- 
ment of  the  object.  A simple  explanation  of  the  technique  follows: 


A hologram,  which  Is  a photograph  containing  Information  about  In- 
tensity and  phase  of  light  reflected  by  an  object,  Is  made  when  a high 
resolution  photographic  plate  Is  exposed  to  two  beams  of  coherent  light 
waves  Intersecting  at  the  surface  of  the  photographic  plate.  One  of  the 
beams  Is  classified  as  the  reference  beam  and  comes  unobstructed  from  the 
laser.  The  other  beam,  known  as  the  object  beam,  comes  from  the  laser  but 
Is  first  reflected  off  an  object  before  reaching  the  photographic  plate. 
The  fact  that  the  light  waves  are  coherent  (light  waves  travel- 
ing through  space  are  of  the  same  frequency  and  phase)  enables  a grating 
to  be  formed  on  the  plate.  This  grating  Is  formed  from  the  constructive 
and  destructive  Interference  of  the  two  light  beams,  l.e. , object  beam 
and  reference  beam.  After  the  exposed  plate  has  been  developed,  returned 
to  Its  original  position,  and  Illuminated  by  a coherent  beam,  the  grating 
on  the  plate  will  cause  the  light  to  be  bent  or  diffracted  such  that  the 
object  beam  Is  recreated.  This  beam,  which  was  originally  diffracted  off 
the  object,  reconstructs  the  Image  of  the  object  appearing  In  Its  original 
position  In  space,  provided  all  of  the  factors  are  the  same  as  In  the 
original  exposure.  In  order  to  view  the  reconstructed  Image,  one  looks 
through  the  hologram  and  sees  the  image  appearing  In  space.  This  Is 
known  as  transmission  holography. 

In  using  the  technique  for  vibration  analysis,  a transmission  holo- 
gram Is  made  of  the  object  to  be  analyzed  (a  disk  In  this  case)  while  It 
Is  stationary.  The  hologram  Is  developed  and  placed  In  Its  original 
position  and  one  observes,  while  looking  through  the  plate,  the  original 
object  and,  superimposed  on  It,  the  Image  of  the  object  reconstructed  by 
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the  hologram.  If  the  object  Is  vibrating,  the  resulting  displacement 
will  cause  fringes  to  appear.  These  are  caused  by  the  light  rays  from 
the  object  Interfering  with  those  from  the  Image,  and  correspond  to  the 
displacement  of  the  object  from  Its  original  stationary  position.  When 
the  object  undergoes  a sinusoidal  forced  vibration,  fringes  will  be  seen 
on  the  object  only  when  all  points  on  the  object  are  vibrating  In  phase 
or  180°  out  of  phase  with  each  other.  Such  Is  the  case  when  an  elastic 
body  Is  vibrating  at  one  of  Its  natural  frequencies.  Thus,  by  sweeping 
through  the  frequency  range  of  Interest,  the  natural  frequencies  and  mode 
shapes  can  be  determined  by  observing  the  frequencies  at  which  fringes 
occur  and  how  they  are  distributed.  This  type  of  testing  Is  termed  "real- 
time"  holography,  In  that  the  fringes  appear  or  disappear  as  a result  of 
the  sinusoidal  forced  vibration. 

Once  the  natural  frequencies  of  the  object  are  determined  using 
"real-time"  holography,  permanent  photographs  of  the  mode  shapes  can  be 
recorded  using  time-averaged  holography.  This  method  consists  of  expos- 
ing the  holographic  plate  with  reflected  coherent  light  from  the  object 
while  It  Is  vibrating  at  one  of  Its  natural  frequencies.  The  resulting 
hologram,  after  proper  development,  shows  the  object  to  be  covered  with 
Interference  fringes.  These  fringes  are  caused  by  Interference  when  the 
object  Is  at  Its  two  extreme  oscillatory  excursions.  The  points  on  the 
object  which  have  no  relative  motion,  l.e. , the  nodal  lines,  reflect  the 
most  light  and  appear  as  the  brightest  fringes  In  the  hologram.  As  we 
travel  outward  from  the  nodal  lines,  each  fringe  encountered  represents 
a contour  of  constant  displacement.  Hence,  the  time-averaged  hologram 
yields  a contour  map  of  the  normal  displacement.  This  can  be  easily  pho- 
tographed by  using  common  high  contrast  photographic  films. 

The  experimental  set-up  used  to  determine  the  natural  modes  and 
frequencies  of  the  steel  disk  are  shown  In  Figures  2 and  3.  A HeNe  50 
MW  laser  (Spectra  Physics  125A)  was  used  as  a light  source.  As  shown  In 
Figure  3,  the  laser  beam  Is  reflected  off  mirrors  (M)  and  Is  divided  Into 
two  parts  by  the  beam  splitter  (BS).  The  object  beam  (0)  Is  expanded  by 
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the  beam  expander  (BE).  The  beam  then  passes  to  another  beam  splitter 
which  reflects  the  light  to  the  disk.  The  purpose  of  this  beam  splitter 
Is  to  make  the  light  waves  perpendicular  to  the  disk  so  that  optimum  Illum- 
ination, especially  under  rotating  conditions,  will  be  achieved  by  using 
retro-reflective  paint  applied  to  the  surface  of  the  disk.  The  object 
beam  Is  then  reflected  from  the  disk  and  onto  the  photographic  plate  (Agfa 
Gevaert  10E75).  Finally,  the  reference  beam  (R)  via  the  beam  splitter, 
mirrors,  and  expanders  goes  directly  to  the  photographic  plate.  The  Inter- 
ference of  the  two  beams  at  the  photographic  plate  enables  a hologram  (real 
time  or  time-averaged)  to  be  made  of  the  circular  disk. 

3.  EXPERIMENTAL  PROCEDURES  - ROTATING  CONDITIONS 

After  determining  the  mode  shapes  and  frequencies  under  static 
conditions,  the  disk  was  mounted  to  a 3/4  HP  tool  post-grinder  with  a 
speed  capability  ranging  from  zero  to  10,000  rpm.  The  first  five  dia- 
metral modes  of  vibration  (In,  Os)  through  (5n,  Os)  and  the  second  couple 
mode  of  vibration  (2n,  Is)  were  determined  using  image  derotated  double- 
pulsed  holography  as  a function  of  speed.  This  method  Incorporates  an 
Image  derotator,  developed  by  P.  Waddell  (Reference  6),  to  compensate 
optically  for  the  circular  motion  of  the  object  giving  It  the  appearance 
of  a stationary  object.  Stetson  (Reference  3)  applied  the  Image  derotator 
In  conjunction  with  standard  double-pulsed  laser  holography  and  demon- 
strated that  Interferometric  holograms  could  be  recorded  of  objects  rotat- 
ing at  high  speeds.  (This  work  was  funded  by  the  Air  Force  Aero-Propulsion 
Laboratory  under  Contract  F3361 5-75-C-201 3 . ) The  work  culminated  In  the 
delivery  and  Installation  of  an  Image  derotator  system  at  Wrlght-Patterson 
Air  Force  Base.  It  was  this  Instrument  which  was  used  to  record  the  vibra- 
tion characteristics  of  the  spinning  disk.  An  Important  advantage  of  this 
system  over  stroboscopic  methods  for  stopping  the  rotation  of  an  object 
Is  Its  ability  to  provide  continuous  motion  compensation  of  the  object. 

This  Is  Independent  of  the  size  of  the  object  and  Is  limited  only  by  the 
mechanics  of  rotating  the  prism,  so  that  tolerances  on  laser  pulse  width 
and  pulse  repetition  rate  are  greatly  relaxed. 
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As  mentioned  previously,  a double-pulsed  ruby  laser  (Apollo  Model 
22HD)  was  used  In  conjunction  with  the  Image  derotator  to  obtain  holo- 
grams under  conditions  of  rotation.  The  fringes  formed  represent  contours 
of  constant  vibratory  displacement.  It  should  be  mentioned  that  P.  Waddell 
(Reference  2)  has  made  attempts  to  use  a brute  force  method  of  gating  a 
laser  Into  one  pulse  per  revolution  so  that  holographic  interferometry 
could  be  accomplished  without  an  Image  derotator.  Extreme  difficulty  Is 
encountered  In  maintaining  registration  of  the  hologram  and  object  (l.e., 
shaft  whirl).  This  together  with  the  low  laser  energy  Inherent  to  the 
short  pulses  required  to  stop  the  rotary  motion,  make  this  method  extreme- 
ly unattractive.  However,  by  Incorporating  an  Image  derotator,  the  object 
appears  stationary  and  the  problems  of  maintaining  registration  of  the 
holograms  over  a cycle  of  revolution  and  the  problems  associated  with  gat- 
ing the  laser  are  eliminated. 


The  experimental  set-up  used  to  determine  the  mode  shapes  and 
frequencies  under  rotating  conditions  Is  shown  In  Figure  4.  This  Is 
a standard  off-axis  holographic  set-up.  A Q-swItched  pulsed  ruby  laser 
Is  used  for  generating  the  laser  light.  It  passes  through  a beam  ex- 
pander and  onto  a beam  splitter.  The  light  that  Is  transmitted  becomes 
the  object  beam,  reflecting  through  a large  double  folding  mirror  assem- 
bly (Ml)  which  directs  It  to  the  disk.  By  virtue  of  retro-reflective 
paint  on  the  disk,  the  beam  Is  reflected  back  and  retraces  the  same  path 
back  to  the  beam  splitter  (BS)  whe**e  It  Is  directed  through  the  derota- 
tor onto  the  holographic  film.  The  reflected  light  from  the  beam  split- 
ter acts  as  the  reference  beam.  It  Is  reflected  off  mirror  M2  onto 
mirror  M3  and  then  onto  the  holographic  film. 

The  Image  derotator  shown  In  Figure  5,  which  Is  the  most  complex 
portion  of  the  holographic  set-up.  Is  comprised  of  a transmissive  folded 
Abbe  prism  (Supplied  by  Optical  Systems  and  Technology,  Inc.,  East  Hart- 
ford, Connecticut)  which  Is  located  within  a hollow  shaft  motor  comprised 
of  an  air  bearing  and  a brushless  DC  torque  motor.  A code  wheel  that  has 
a number  of  equally  spaced  circumferential  markers  Is  attached  to  the 


derotator  shaft.  A light-emitting  diode  (LED)  Is  used  to  detect  the 
markers  on  the  code  wheel  and  thus  monitor  the  rotational  speed  and  phase 
of  the  motor.  Another  code  wheel  shown  In  Figure  6 Is  fixed  to  the  shaft 
of  the  lathe  motor  which  drives  the  disk.  It  also  has  a LED  pickup.  The 
signal  which  is  received  from  this  pickup  Is  compared  to  the  derotator 
signal  via  a control  box  unit.  This  unit  adjusts  the  voltage  to  the  de- 
rotator motor  so  that  Its  rotational  rate  Is  half  that  of  the  rotating 
disk  and  Is  phase-locked  with  It.  Also  shown  Is  an  electromagnet  used 
to  excite  the  disk  and  a capacitance  pickup  used  to  detect  resonance. 


Because  the  derotator  Is  traveling  at  half  the  angular  velocity  of 
the  object,  continuous  motion  compensation  Is  achieved.  An  explanation 
of  the  principles  Involved  follows.  Figure  7 Is  a diagram  of  a folded 
Abbe  prism,  which  Is  comprised  of  three  prisms  mated  together,  and  Is 
located  in  the  center  of  the  derotator.  Note,  the  top  and  bottom  surfaces 
of  the  prism  system  are  mirrored  and  thus  reflect  the  light  waves.  Also, 
depending  on  the  relative  angle  In  which  the  light  waves  strike  the  Inter- 
facing surfaces  of  the  prisms,  the  light  Is  either  reflected  or  trans- 
mitted. Now,  If  an  object  (note  arrow  In  vertical  position)  Is  viewed 
through  the  wedged  prism  system.  It  appears  to  be  Inverted.  Rotate  the 
object  180°  in  the  clockwise  direction  and  the  derotator  90*  In  the  same 
direction,  remembering  the  derotator  travels  at  half  the  speed  of  the 
object,  the  object  as  viewed  through  the  prism  remains  stationary.  Again, 
rotate  the  object  another  180*  and  correspondingly  the  derotator  90®. 
Observe  the  arrowhead  remains  In  the  same  position.  Therefore,  for  any 
angular  position  of  the  object.  It  appears  stationary  as  viewed  through 
the  derotator  prism  system,  provided  the  prism  system  Is  rotated  at  half 
the  angular  speed  of  the  object.  This  Is  the  principle  on  which  the 
derotator  works,  A more  thorough  explanation  of  the  system  Is  found  In 
Reference  3. 
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The  disk  shown  In  Figures  7 and  8 was  driven  by  means  of  a 3/4  HP 
tool  post-grinder  motor  connected  by  a belt  to  a two-bearing  Internal 
grinding  spindle.  The  disk  was  machined  with  a 0.50  Inch  diameter  hole 
for  sliding  onto  a half-inch  drive  spindle.  A keyway  was  used  to  prevent 
slippage  of  the  disk  relative  to  the  spindle.  A 3/8  Inch-16  nut  applied 
with  a torque  of  100  Inch-lb  fixed  the  disk  to  the  spindle.  Fixed  to  the 
opposite  end  of  the  spindle  was  the  code  wheel.  This  drive  system  had  a 
capability  of  10,000  rpm  which  Is  the  upper  limitation  of  the  derotator. 
Prior  to  spinning  the  disk,  Its  burst  speed  was  calculated  to  ensure 
operator  safety.  The  results  (see  Appendix  A)  Indicates  a safety  margin 
of  approximately  two  which  was  considered  sufficient.  The  drive  system 
was  more  than  adequate  to  achieve  a true  center  of  rotation  of  the  disk, 
a requirement  necessary  to  provide  continuous  alignment  between  the  disk 
axis  of  rotation  and  the  derotator  optical  axis  of  rotation.  The  rota- 
tional speed  of  the  grinder  was  controlled  by  a simple  rheostat. 


To  excite  the  disk  under  conditions  of  rotation  Into  one  of  Its 
resonant  modes,  either  an  electromagnet  or  siren  was  used.  It  was  placed 
a fixed  distance  parallel  to  the  disk  surface  at  Its  outer  circumference 
and  driven  by  a sweep  oscillator.  A stationary  vibration  pickup,  typi- 
cally a capacitance  probe,  was  placed  near  the  back  face  of  the  disk  to 
monitor  out-of-plane  displacement.  As  In  the  nonrotating  configuration, 
the  excitation  frequency  signal  and  the  frequency  signal  from  the  capa- 
citance pickup  were  fed  respectively  to  the  horizontal  and  vertical  plates 
of  an  oscilloscope.  The  resulting  Llssajous  patterns  were  used  to  deter- 
mine resonance.  Figure  9 shows  the  Instrumentation  used  to  establish 
this  condition.  Holographic  Interferograms  using  the  derotator  In  con- 
junction with  the  ruby  laser  were  made  at  various  resonant  conditions  as 
a function  of  speed.  The  excitation  frequency  and  disk  rpm  were  recorded 
for  each  of  these  Interferograms.  In  this  manner,  a frequency- speed  dia- 
gram showing  the  various  resonant  conditions  for  the  spinning  disk  was 
generated. 
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Two  factors  should  be  noted  In  utilizing  the  Image  derotator  for 
execution  of  holographic  Interferometry  as  described  above.  First,  the 
reference  beam  does  not  pass  through  the  derotator  and  hence  remains 
linearly  polarized  throughout  Its  path  of  travel.  The  object  beam  by 
virtue  of  passing  through  the  derotator  undergoes  a change  In  polariza- 
tion. To  record  a hologram,  this  light  must  have  the  same  polarization 
as  the  intended  reference  beam.  Fortunately,  as  noted  by  Stetson  (Refer- 
ence 3),  the  construction  of  the  folded  Abbe  Inverting  prism  afforded  the 
characteristics  of  a half-wave  retarding  plate.  By  simply  placing  a quar- 
ter-wave retarder  plate  at  the  derotator  entrance,  the  reflected  linearly 
polarized  object  beam  Is  converted  to  circular  polarized  light.  The  beam 
is  converted  back  to  linearly  polarized  light  with  another  quarter-wave 
retarder  after  the  derotator.  Thus,  It  Is  possible  to  achieve  nearly 
complete  polarization  compensation. 


Secondly,  since  the  object  is  Illuminated  with  linearly  polarized 
light,  under  conditions  of  rotation,  the  reflected  light  will  undergo  a 
cyclic  change  In  orientation  relative  to  the  derotated  Image  of  the  disk. 
In  other  words,  the  Image  appears  stationary  but  the  linearly  polarized 
light  Is  rotating  at  the  angular  velocity  of  the  disk  and  In  the  opposite 
direction.  There  would  be  a total  lack  of  correlation  and  hence  no  Inter- 
ferogram  if  the  disk  were  allowed  to  rotate  more  than  90°  between  pulses 
of  the  ruby  laser.  This  was  easily  accomplished  by  limiting  the  laser 
pulse  width  to  less  than  50  us.  At  10,000  rpm,  this  would  amount  to 
less  than  3°  of  disk  rotation. 

4.  ALIGNMENT 

Prior  to  recording  double  exposure  holograms  of  the  disk  resonant 
response  under  conditions  of  rotation,  optical -mechanical  alignment  of 
the  system  has  to  be  achieved.  Initially,  a point  light  source  from  a 
HeNe  laser  Is  directed  onto  the  object  via  the  same  optics  used  with 
the  ruby  laser.  The  derotator  Is  switched  to  manual  and  set  to  rotate 
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at  a slow  speed.  If  the  axis  of  rotation  of  the  derotator  Is  pointed 
directly  at  the  point  source,  the  Image  formed  will  describe  a circle. 

If  the  derotator  prism  axis  Is  not  colllnear  with  the  rotational  axis 
of  the  drive  motor,  the  Image  formed  will  describe  an  epicycle.  The 
epicycle  Is  reduced  to  a circle  by  adjusting  the  angular  orientation  of 
the  Abbe  prism.  Now,  It  Is  necessary  to  reduce  the  circle  to  a single 
point.  This  will  ensure  that  the  derotator  optical  axis  and  the  object 
rotational  axis  are  now  colllnear  and  aligned.  This  Is  accomplished  by 
adjusting  the  mirror  system  Ml  In  Figure  4 with  the  object  stationary 
and  then  a final  adjustment  with  it  spinning  and  phase-locked  with  the 
derotator.  A check  of  the  alignment  Is  made  under  rotating  conditions 
by  making  a double  exposure  hologram  without  any  resonant  excitation. 

The  resulting  exposed  hologram  will  yield  a group  of  parallel  bias  fringes 
across  the  face  of  the  disk  surface  whose  spacing  is  proportional  to  the 
degree  of  misalignment.  If  the  number  of  these  bias  fringes  approaches 
the  same  order  of  fringes  formed  from  disk  normal  displacement  resulting 
from  resonant  excitation,  a confusing  hologram  will  result  and  the  sys- 
tem has  to  be  realigned.  However,  If  the  double  exposure  hologram  has 
fewer  than  five  bias  fringes,  the  system  degree  of  alignment  Is  more  than 
adequate  to  achieve  distinctly  defined  resonant  mode  shapes  under  rotating 
conditions. 

5.  NUMERICAL  ANALYSIS 

The  numerical  portion  of  the  study  was  based  on  finite  element 
methods  and  classical  solutions  readily  available  for  circular  plates. 

Only  the  first  five  diametral  modes  and  the  coupled  (2n-ls)  mode  were 
Investigated. 

Lelssa  (Reference  7)  has  generated  tables  of  eigenvalues  for  annu- 
lar plates  consisting  of  circular  outer  boundary  (a)  and  a concentric 
circular  Inner  boundary  (b)  with  constant  thickness  (T).  Various  com- 
binations of  simple  boundary  conditions  are  presented  for  a Poisson's 
ratio  of  1/3.  These  tables  were  applied  to  a disk  which  was  assumed  to 
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be  free  on  the  outside  and  clamped  on  the  inside.  The  outside  radius 
was  5 Inches  and  Inside  radius  one  half  inch  with  a thickness  of  1/16 
inch.  This  simulates  the  disk  geometry  which  was  used  In  the  experi- 
mental set-up.  Based  on  this  configuration  and  Lelssa's  tables,  the 
first  five  nodal  diameters  and  the  coupled  (2n-ls)  mode  were  calculated 
for  nonrotating  conditions. 

The  finite  element  computer  program  NASTRAN  was  used 
to  compute  the  natural  frequencies  and  mode  shapes  for  a circular  annu- 
lar plate  under  nonrotating  conditions  and  Increasing  centrifugal  loading 
up  to  10,000  rpm.  The  annulus  of  the  plate  was  fixed  against  translation 
and  rotation  with  the  geometry  being  the  same  as  that  used  in  the  classi- 
cal solutions  and  experimental  set-up. 

The  type  of  element  used  was  the  quadrilateral  plate  bending 
element  called  CQUAD2  (this  is  a name  for  a 4-sided  element  in  the  finite 
element  method  used  in  NASTRAN).  This  element  has  both  inplane  and  bend- 
ing stiffness  that  assumes  a solid  homogeneous  cross  section.  The  ele- 
ments used  were  kept  as  nearly  square  as  practicable  because  the  accuracy 
tends  to  deteriorate  as  the  aspect  ratio  of  the  quadrilateral  Increases. 

The  mesh  composed  of  CQUAD2  elements  was  used  in  the  eigenvalue 
extraction  routine  of  the  Inverse  power  method  to  determine  the  natural 
frequencies  of  vibration.  This  Is  a stand-alone  method  In  NASTRAN  and 
the  routine  finds  all  the  eigenvalues  within  a domain  specified  by  the 
user. 

In  addition  to  computing  natural  frequencies,  the  version  of  NASTRAN 
used  at  the  Air  Force  Aero-Propulsion  Laboratory  permits  contour  plots 
to  be  made  of  the  vibration  mode  shapes  showing  out-of-plane  displace- 
ment in  the  Z-di recti on  (normal  to  the  plane  of  the  plate).  This  plot- 
ting was  conducted  for  the  modes  of  Interest. 
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In  this  study  we  wanted  to  investigate  the  vibratory  response  of 
the  plate  under  the  additional  Influence  of  centrifugal  loading.  This 
was  necessary  to  verify  the  results  obtained  experimentally  using  the 
Image  derotator.  Because  of  the  cost  of  running  such  an  analysis  on  the 
computer,  only  one  mode  of  vibration  was  Investigated  at  rpms  of  2,000, 
4,000,  6,000,  and  8,000. 

Rigid  format  13  of  NASTRAN  was  used  to  compute  the  natural  fre- 
quencies and  mode  shapes  for  the  one  case.  NASTRAN  performs  the  vibra- 
tional analysis  in  a two-step  operation.  First,  it  takes  the  value  of 
the  angular  velocity  Input  to  it  and  creates  a set  of  forces  that  act 
on  the  element  centroids  of  the  finite  element  model.  These  forces  are 
a function  of  the  mass  of  the  element,  their  locations  with  respect  to 
the  axis  of  rotation,  and  the  angular  velocity.  Using  these  forces,  a 
static  analysis  Is  carried  out  (rigid  format  1)  and  the  resulting  stresses 
are  used  to  generate  a centrifugal  stiffness  matrix.  In  the  second  step 
of  this  operation,  the  centrifugal  stiffness  matrix  Is  added  to  the  NAS- 
TRAN-g enerated  elastic  stiffness  matrix  and  an  eigenvalue  analysis  Is 
carried  out  (rigid  format  3)  yielding  the  natural  frequencies  of  the 
structure.  This  two-step  process  Is  carried  out  automatically  using 
rigid  format  13  of  NASTRAN. 
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SECTION  III 
RESULTS 

1.  EXPERIMENTAL  RESULTS 

In  order  to  have  a reference  for  the  disk  mode  shapes  obtained 
under  rotating  conditions,  the  first  five  diametrical  modes  of  vibration 
(In,  Os)  through  (5n,  Os)  and  second  coupled  mode  of  vibration  (2n,  Is) 
for  the  steel  disk  fixed  at  its  center  were  determined  using  time-averaged 
holography.  Except  when  noted,  the  Os  portion  of  the  mode  notation  shall 
be  Implied. 

Figures  10  and  11  show  the  six  resonant  modes  obtained  statically 
using  siren  excitation  and  a standard  off-axis  holographic  set-up  with  a 
50  MW  HeNe  laser  as  described  previously.  Resonance  was  determined  by 
monitoring  the  out-of-plane  displacement  with  a capacitance  probe  while 
varying  the  frequency  via  a sweep  oscillator.  The  outputs  of  each  sig- 
nal were  fed  to  the  horizontal,  and  vertical  Inputs  to  an  oscilloscope, 
with  resonance  being  determined  by  the  resulting  Lissajous  patterns. 

The  brightest  fringes  in  the  figures  for  the  holograms  are  the 
node  lines  which  are  points  of  zero  displacement  forming  diametrically 
across  the  face  of  the  disk;  hence,  the  terminology  diametral  modes. 

The  couple  mode  shown  In  Figure  11  (1052  HZ)  is  a combination  of  dia- 
metrical modes  plus  circumferential  modes  running  around  the  periphery 
of  the  disk. 


The  node  line  orientation  was  governed  primarily  by  the  location 
of  the  siren  relative  to  the  disk.  The  point  of  excitement  could  be 
localized  to  approximately  a one-inch  diameter  circle.  At  this  point, 
the  mode  pattern  oriented  Itself  such  that  an  antinode  (point  of  maxi- 
mum displacement  which  bisects  two  adjacent  node  lines)  was  yenerated 
at  the  siren  excitation  location.  However,  It  was  discovered  that  by 
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lowering  the  siren  excitation  power  the  modes  assumed  a preferred  orien- 
tation. This  was  comprised  of  two  sets  of  node  patterns  for  each  mode 
of  vibration.  They  were  essentially  the  Inverse  of  each  other  with  the 
nodal  lines  of  one  becoming  the  anti  nodal  lines  of  the  other;  l.e.,  the 
node  line  which  was  the  point  of  zero  displacement  now  becomes  the  anti- 
node  line  or  point  of  maximum  displacement.  The  preferred  orientation 
Is  a result  of  the  imperfection  of  the  disk  (Reference  4),  which  is  pres- 
ent in  all  practical  cases.  The  two-mode  patterns  differed  In  natural 
frequencies  by  an  amount  depending  on  the  degree  of  imperfection.  These 
modes  are  quite  independent  of  each  other,  and  their  nodal  patterns  coin- 
cide with  those  of  maximum  and  minimum  frequency.  The  nodal  diameter  of 
one  and  an  antlnodal  diameter  of  the  other  will  pass  through  the  Imper- 
fection. The  frequency  of  the  former  will  be  greater  than  that  of  the 
latter. 

For  the  steel  disk  In  this  experiment,  the  frequency  spread  ranged 
from  0.5  to  1.5  hertz  depending  on  the  mode  being  excited.  This  small 
difference  Is  indicative  of  the  high  degree  of  radial  symmetry  present 
in  the  mass  and  stiffness  characteristics  of  the  disk.  It  can  be  shown 
that  In  the  case  of  the  disk  with  perfect  radial  symmetry,  that  each  of 
the  mode  pairs  will  degenerate  to  a single  mode  shape  and  corresponding 
resonant  frequency  having  random  node  orientation.  Thus,  a perfect  disk 
would  have  no  preferred  node  orientation.  An  Imperfect  disk  under  condi- 
tions of  free  vibration  would  have  two  node  patterns  for  each  mode  whose 
frequency  spread  depends  on  the  amount  of  unsymmetry.  Under  conditions 
of  forced  vibration,  only  one  mode  pattern  will  appear  with  Its  orien- 
tation depending  on  the  excitation  source  location. 

The  disk  was  mounted  to  the  lathe  motor  and  Its  resonance  response 
was  determined  under  rotating  conditions.  Figures  12  and  13  show  the 
results  which  describe  the  phenomena  mentioned  previously  In  the  discus- 


sion. This  demonstrates  that  for  each  mode  of  vibration  two  excitation 
frequencies  will  excite  the  same  resonant  mode.  In  other  words,  at  a 
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specific  disk  rpm,  the  2n  mode  of  vibration  will  be  excited  into  reso- 
nance by  two  different  excitation  frequencies,  one  that  is  2ft  greater 
than  the  2n  natural  frequency  and  one  that  Is  2ft  less  than  the  2n  natur- 
al frequency.  Shown  are  the  upper  and  lower  resonant  conditions  for  the 
In,  2n,  3n,  4n,  and  5n  modes,  respectively. 

The  plots  (Figures  12  and  13)  depict  the  frequency  of  an  electro- 
magnetic exciter  required  to  generate  a resonant  disk  mode  versus  the 
disk  speed.  The  data  points  represent  a maximum  response  indication 
from  a capacitance  probe  placed  near  the  back  face  of  the  disk  (see 
Figure  4,  "Vibration  Pickup").  The  curve  joining  each  set  of  data  points 
is  a second  order  least  squares  fit.  As  discussed,  the  upper  curve  for 
each  mode  of  vibration  is  a plot  of  pp  + n ft.  As  the  disk  speed  decreases 
to  zero,  the  upper  and  lower  branches  of  the  excitation  frequency  con- 
verge to  the  corresponding  resonant  frequency  for  the  stationary  disk. 

The  dashed  lines  in  Figure  12  represent  extrapolated  data  and  show  the 
2n  and  3n  modes  crossing  the  abscissa  at  5600  and  7600  rpm,  respectively. 

At  these  points,  a resonant  mode  can  be  excited  by  a static  force.  This 
would  be  a particularly  dangerous  region  of  operation  for  a turbine  engine 
disk  since  a circumferential  variation  In  the  existing  pressure  field  could 
generate  a resonant  response.  These  two  resonant  points  are  often  re- 
ferred to  In  a "Campbell  Diagram"  used  by  turbine  engine  designers.  The 
diagram  is  established  by  bisecting  the  upper  and  lower  resonant  lines 
shown  In  Figures  12  and  13.  This  determines  the  natural  frequencies  for 
each  disk  mode.  Then  engine  order  lines  are  plotted  which  are  equal  to 
the  number  of  excitations  per  disk  revolution  (IE,  2E,  3E,  etc.).  Now, 
the  2E  engine  order  line  will  cross  the  2n  mode  natural  frequency  at  5600 
rpm  and  have  a value  of  5600  x 2E/60  = 187  Hz.  Likewise,  the  3E  engine 
order  line  will  cross  the  3n  mode  natural  frequency  at  7600  rpm  and  have 
a value  of  7600  x 3E/60  = 380  Hz.  The  Campbell  Diagram  Is  then  used  by 
the  designer  to  Identify  possible  resonant  points  which  should  be  avoided 
during  engine  operation. 
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Figures  14  and  15  show  the  same  disk  resonant  mode  shapes  that 
were  obtained  under  stationary  conditions  (Figures  10  and  11)  obtained 
using  Image  derotated  holographic  interferometry  at  specified  disk  speeds. 
The  mode  shapes  shown  In  the  figures  are  double  exposure  holograms  made 
using  the  pulsed  ruby  laser  with  a pulse  separation  of  20  us.  It  Is 
seen  that  the  fringe  contrast  is  quite  good  and  almost  entirely  propor- 
tional to  disk  out-of-plane  displacement  with  very  little  distortion 
Introduced  from  bias  fringes  caused  by  opto-mechanlcal  misalignment. 

Each  mode  shape  In  Figures  14  and  15  corresponds  to  some  point  on  one 
of  the  excitation  frequency  - disk  speed  curves  in  Figures  12  and  13. 

For  example,  the  4n  mode  shape  In  Figure  15  was  excited  at  724  Hz  when 
the  disk  was  rotating  at  2620  rpm.  This  point  lies  on  the  upper  branch 
of  the  4n  curve  In  Figure  13.  The  actual  natural  frequency  for  the  mode 
at  this  rpm  would  thus  be: 

P4  = Fe4  - n n = 724  - (4)  (2620)/60  = 549  Hz 

This  Is  an  Increase  of  12  Hz  over  the  4n  natural  frequency  for  the 
nonrotating  disk.  Figure  16  shows  the  same  4n  mode  shape  obtained  while 
the  disk  was  running  at  7920  rpm.  The  very  minor  contribution  of  bias 
fringes,  l.e.,  the  two  fringes  running  approximately  across  the  diameter 
of  the  disk,  demonstrates  the  stability  of  the  system  even  at  speeds 
near  the  upper  design  limit  of  the  derotator  (10,000  rpm). 


Figures  14,  15,  and  16  show  modes  of  vibration  that  are  fairly 
pure  In  their  nature,  l.e.,  mode  shapes  characteristic  of  classical 
disk  vibration  theory.  Disk  mode  shapes  were  also  recorded  that  were 
apparently  combinations  of  more  than  one  resonant  mode  shape.  This 
type  of  response  would  manifest  Itself  as  the  actual  resonant  mode  shape 
being  excited  and  distorted  typically  by  the  In  mode  shape.  This  type 
of  response  was  thought  to  be  caused  by  an  additional  disk  excitation 
force  resulting  from  rotor  dynamic  imbalance  of  the  disk  shaft.  This 
phenomenon  was  subsequently  verified  by  carrying  out  a power  spectral 
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density  analysis  of  the  vibration  transducer  signal  with  and  without 
disk  excitation.  The  excitation  component  due  to  shaft  Imbalance  was 
not  found  to  present  a problem  provided  that  the  amplitude  of  the  elec- 
tromagnetic excitation  was  an  order  of  magnitude  larger  than  the  force 
from  the  dynamic  shaft  Imbalance.  There  was  no  added  complication  In 
Increasing  the  electromagnetic  force  amplitude  since  the  laser  pulse 
separation  could  be  decreased  proportionally  In  order  to  avoid  "wash  out" 
of  the  holographic  fringes  from  too  large  a disk  displacement. 

2.  FINITE  ELEMENT  ANALYSIS  AND  THEORETICAL  CORRELATIONS 

The  finite  element  results  were  based  on  a mesh  shown  In  Figure 
17  which  was  generated  to  simulate  the  steel  disk  used  In  the  experimental 
work.  It  was  internally  clamped  along  an  inner  radius  (b  = 1/2  Inch) 
and  unconstrained  along  the  outer  radius  (a  = 5 Inches).  Care  was  taken 
to  maintain  a length  to  width  ratio  of  approximately  one  for  the  element 
blocks  to  ensure  accuracy. 

The  finite  element  results  are  shown  in  Figures  18  and  19  for  the 
modes  of  vibration  which  were  determined  experimentally  under  static 
conditions.  The  contour  numbers  In  the  figures  correspond  to  values  of 
displacement  in  the  Z-direction  (out  of  the  paper)  such  that  the  integers 
1,  2,  3,  4,  5,  6,  7,  8,  9,  10,  and  11  refer  to  the  Z-dlsplacement  values 
-1.,  -0.8,  -0.6,  -0.4,  -0.2,  0,  1.3,  0.4,  0.6,  0.8,  and  1.,  respectively. 
Note  that  contour  number  6 corresponds  to  a Z-dlsplacement  of  zero  and 
hence,  defines  the  vibration  nodal  lines  for  each  mode  shape.  For  the 
purpose  of  this  work,  the  nodal  line  6 was  further  darkened  by  hand  to 
enhance  Its  visibility.  The  contour  plots  shown  serve  to  graphically 
display  the  deformation  characteristics  of  each  mode  of  vibration.  Note 
that  all  the  modes  exhibit  either  diametral  and/or  circumferential  symme- 
try about  the  surface  of  the  disk  which  Is  typical  of  classical  plate 
vibration  response. 
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Lelssa's  tables  based  on  classical  solutions  for  circular  plates 
were  used  to  establish  the  nondlmens Iona  1 frequency  parameter  X2  where 

X2  = (joa2(p/D)1/2 

(4) 

0 = ET3/ 1 2 (1-v2) 

2 

Using  the  tabulated  values  of  X , we  solved  Equation  4 for  to  In  radians/ 
sec  for  the  same  frequencies  determined  experimentally.  For  this  case, 
the  following  properties  (same  as  the  steel  disk)  were  used:  radius,  a = 
5.0  Inches,  thickness,  T = 0.0625  Inches,  Poisson's  ratio,  v = 0.3, 

Young's  modulus,  E = 3.05  x 107  psl,  mass/area,  p=  1.768  x 10“2  lbm/1n2. 
The  parameter  D Is  known  as  the  flexural  rigidity  and  is  constant  for  a 
given  configuration.  Table  1 shows  the  correlation  between  the  experi- 
mental results,  finite  element,  and  frequencies  determined  using  Lelssa's 
tables  (Reference  7).  There  was  a large  discrepancy  between  the  numerical 
techniques  and  the  experimental  results.  The  difference  was  judged  to  be 
caused  primarily  by  the  boundary  conditions  imposed  on  the  numerical 
technique  versus  the  steel  hub  configuration  used  In  the  experimental 
work.  An  experiment  was  conducted  (see  Appendix  B)  using  a different 
configuration  for  clamping  the  disk  which  was  considered  to  be  more  repre- 
sentative of  the  numerical  boundary  constraints.  The  disk  was  tested  ex- 
perimentally and  results  compared  to  analytical  calculations  of  frequen- 
cies. Correlation,  in  most  cases,  was  less  than  a 21  difference,  leading 
to  the  conclusion  that  the  Integral  hub  for  the  steel  disk  does  not  con- 
strain the  Inner  radius  as  was  originally  thought.  However,  the  purpose 
of  the  research  was  to  prove  the  merits  of  the  Image  derotator,  so  that 
discrepancies  between  experimental  and  analytical  techniques  are  of  little 
concern. 


The  finite  element  program  was  further  used  to  provide  an  analyti- 
cal check  on  the  effect  of  centrifugal  loading  on  the  resonant  frequencies 
of  the  disk.  Figure  20  shows  a plot  of  the  disk  mode  (3n,  Os)  resonant 
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frequency  as  a function  of  rpm  for  both  the  finite  element  and  experi- 
mental results.  The  experimental  resonant  frequencies  were  approximately 
4%  below  those  computed  from  the  finite  element  model.  This  was  felt  to 
be  sufficiently  accurate  considering  the  large  effects  the  hub  can  have 
on  the  boundary  conditions  as  previously  determined  for  the  nonrotating 
calculations. 
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SECTION  IV 
CONCLUSIONS 

Image  derotation  techniques  In  conjunction  with  double-pulsed  laser 
technology  has  opened  new  avenues  for  research  In  addressing  the  area  of 
dynamic  structural  response  of  rotating  structures.  The  emphasis  of  the 
research  was  to  demonstrate  the  graphic  results  which  can  be  obtained  by 
applying  image  derotated  holographic  interferometry  to  a simple  disk  un- 
der conditions  of  rotation.  With  this  in  mind,  it  should  be  noted  that 
the  characteristic  movement  of  the  mode  shapes  relative  to  the  disk  re- 
quires a much  more  detailed  analysis  than  was  possible  within  the  scope 
of  this  present  work.  The  shapes  recorded  under  rotating  conditions 
such  as  Figures  14  and  15  do  not  provide  information  as  to  the  relative 
movement  of  the  mode  in  relation  to  the  disk.  One  has  to  determine  whe- 
ther the  mode  is  fixed  to  the  disk  and  hence  rotating  with  it,  or  is 
stationary  in  initial  space  with  the  disk  rotating  through  it,  or  is  a 
movement  of  some  combination  in  between.  It  can  be  easily  shown  that 
the  derotator  and  pulsed  laser  combinations  permit  disk  mode  shapes  to 
be  recorded  under  any  of  these  conditions. 

Attempts  were  made  to  determine  the  relative  movement  of  the  mode 
shapes  by  consecutively  recording  a series  of  interferograms  while  keep- 
ing the  disk  rpm  and  vibration  excitation  amplitude  constant  and  observ- 
ing the  location  of  the  node  lines  relative  to  an  arrowhead  drawn  on  the 
disk.  The  results  were  Inconclusive  because  of  difficulty  encountered 
In  keeping  the  rpm  constant  over  a period  of  time.  This  drifting  caused 
phase  shifts  in  the  vibratory  response  of  the  disk  and  the  resulting 
holograms  showed  a random  orientation  of  the  mode  patterns  and  thus  they 
did  not  seem  fixed  to  the  disk  or  In  Inertial  space.  Tobias  (Reference 
4)  carried  out  a detailed  study  of  the  vibrational  characteristics  of 
an  Imperfect  rotating  disk  and  analytically  demonstrated  that  the  rela- 
tive movement  of  the  mode  shapes  was  directly  related  to  the  phase  shifts 
as  the  rotating  disk  passes  through  a resonant  point.  It  Is  felt  that  If 
the  same  phase  could  be  obtained  over  a given  period  of  time,  the  move- 
ment of  the  node  lines  could  be  determined. 
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The  above  questions  still  need  answers  in  order  to  completely 
understand  the  resonant  response  of  rotating  structures  to  provide  in- 
creased assurance  as  to  their  structural  Integrity.  The  phenomena  of 
vibrational  response  of  a rotating  disk  Is  complex  In  Itself;  the  addi- 
tion of  blades  to  the  structure  further  complicates  the  problem.  This 
is  demonstrated  in  Figure  21,  which  shows  a bladed  disk  assembly  vibrat- 
ing In  one  of  Its  many  resonant  modes  under  nonrotating  conditions  using 
time-averaged  holography.  It  demonstrates  the  complex  mode  shapes  which 
can  be  obtained. 

Through  the  advent  of  image  derotated  holographic  interferometry, 
the  above  questions  will  be  addressed  in  future  research  and  the  complex 
vibratory  behavior  exemplified  by  Figure  21  shall  be  studied. 
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APPENDIX  A 

BURST  SPEED  CALCULATION 

In  the  Interest  of  safety.  It  was  necessary  to  calculate  the  burst 
speed  of  the  steel  disk  which  was  used  for  the  rotational  aspects  of  the 
experimental  Investigation.  Rotational  speeds  approaching  10,000  rpm 
were  expected  and  subsequently  used  as  the  design  limit. 


According  to  Roark  (Reference  8),  formulas  for  calculating  burst 
speeds  should  be  based  on  an  average  stress  instead  of  a maximum  stress. 
This  is  because  the  formulas  for  determining  the  maximum  stress  pre- 
supposes elastic  conditions  and  does  not  account  for  plastic  yielding 
prior  to  failure.  This  yielding  tends  to  equalize  the  stress  intensity 
along  a radial  line. 

Roark's  formula  for  calculating  the  average  stress,  cr  is: 

a 

afl  =0.0008636u)2  (a3  -b3) 

a - b 

where  6 is  the  density,  uj  is  the  angular  velocity,  a and  b,  respectively, 
are  the  outside  and  inside  radii.  Experiments  have  shown  that  for  the 
worst  cases,  annular  disks  have  failed  when  a reaches  61 . 5%  of  the  ulti- 
mate strength.  Based  on  this  criteria,  a value  of  36,900  psl  was  assumed 
for  the  design  stress  and  the  above  equation  solved  for  <d.  This  was 
converted  to  rpm  or  burst  speed  and  found  to  be  23,031  rpm.  This  value 
compared  to  the  design  limit  of  10,000  rpm  for  the  derotator,  was  suf- 
ficient to  provide  a margin  of  safety. 
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APPENDIX  B 

HUB  EFFECTS  ON  FREQUENCIES 

The  discrepancies  found  between  the  numerical  and  experimental 
results  in  determining  modes  of  vibrations  were  largely  attributed  to 
the  integral  hub  machined  for  the  steel  disk.  This  was  verified  by  inves- 
tigating another  means  of  clamping  a circular  disk  to  a shaft  and  again 
comparing  the  experimental  frequencies  with  classical  numerical  solutions. 

A 0.0625-inch-thick  aluminum  disk,  10  inches  in  diameter,  with  a 
0. 25-inch  diameter  center  hole,  was  manufactured  without  an  integral  hub 
as  described  previously  for  the  steel  disk.  A new  clamping  fixture,  as 
shown  in  Figure  21,  was  made  from  aluminum  and  designed  to  simulate  an 
annular  disk  consisting  of  a free  circular  outer  boundary  and  a clamped 
circular  inner  boundary  with  an  inner-to-outer  radius  ratio  of  0.1.  It 
was  felt  that  this  design  would  be  more  representative  than  that  of  the 
Integral  hub  for  matching  the  boundary  constraints  used  in  the  numerical 
calculations,  i.e. , zero  slope  and  zero  deflection  at  the  inner  radius. 

The  aluminum  disk  was  mounted  with  new  clamps  to  an  isolation  table 
and  excited  into  resonance,  using  a siren  driven  by  a sweep  oscillator. 
Real-time  holography  was  used  to  determine  the  first  nine  modes  of  vibra- 
tions which  included  circumferential  modes.  Table  2 is  a listing  of  the 
experimental  results  compared  to  frequencies  determined  analytically  us- 
ing Leissa's  tabulated  values  of  eigenvalues  along  with  the  flexural 
rinidity  constant  determined  for  the  aluminum  disk. 

Excellent  correlation  was  obtained  between  the  experimental  and 
numerical  techniques  for  established  modal  frequencies.  In  most  cases, 
the  difference  was  less  than  2%.  The  larger  variation  found  for  the  (In, 
Os)  mode  can  be  attributed  to  the  fact  that  this  mode  exerts  a couple  on 
the  shaft  transverse  to  its  length  with  the  Inner  radius  not  being  rigidly 
constrained  as  In  the  numerical  solutions.  The  other  modes  are  symmetri- 
cal and  thus  balanced  in  their  reactions  on  the  shaft  resulting  In  zero 
slope  and  zero  deflections  around  the  inner  circumference  of  clamping  as 
in  the  numerical  solution. 
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It  was  concluded  that  the  integral  hub  manufactured  with  the  steel 
disk  has  a large  effect  on  the  stiffening  characteristics  and  boundary 
conditions  of  the  system  which,  in  turn,  affects  the  experimental  results. 
However,  the  purpose  of  the  research  was  to  prove  the  merits  of  image  de- 
rotator In  studying  the  vibrational  characteristics  of  a rotating  struc- 
ture. More  importantly,  one  should  note  that,  even  the  simplest  geometry 
of  a steel  disk  with  an  integral  hub  is  very  difficult  to  model  numerically 
and  results  can  often  be  misleading.  It  is  imperative  that  experimental 
techniques  be  continuously  developed  for  the  Investigations  of  structures, 
especially  more  complex  configurations  such  as  bladed  disk  assemblies. 
Numerical  methods  have  to  have  input  from  experimental  methods  to  ensure 
their  validity. 
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APPENDIX  C 
EQUIPMENT 

The  following  is  a list  of  test  equipment  used  throughout  the 
experimental  aspects  of  the  study: 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 

13. 


Double  Pulse  Ruby  Laser 
Apollo  Lasers,  Model  22,  S/N  283 
PW  20n  Sec,  2.5J  Pulse 

He-Ne  Laser 

Spectra  Physics  Model  125A,  S/N  3109-875 
CW;  50MW  Power 

Storage  Oscilloscope 

Tektronix  Inc.,  Model  7633,  S/N  JH  0 177 
Oscilloscope 

Tektronix  Inc.,  Type  531,  S/N  1051 
Universal  Counter 

Hewlett  Packard,  Model  5325B,  S/N  5090J 
Sweep  Oscillator 

Spectral  Dynamics  Corp,,  Model  SD104A-5,  S/N  1477 

Tracking  Filter  (Dynamic  Analyzer) 

Spectral  Dynamics  Corp.,  Model  SD101B,  S/N  645 

Charge/Voltage  Amplifier 

Unhol tz-Dickie,  Model  D22,  S/N  PJ6501 

Tool  Post  Grinder 

Dumore,  Cat.  No.  57-011,  S/N  8476-01914 
RPM  10,500  HP  3/4 

Power  Amplifier 

Bell  Sound  Division,  Model  Mark  40,  S/N  JC8056 
Amplitude  Servo/Monitor 

Spectral  Dynamics  Corp.,  Model  SD  105C,  S/N  5969J 

35MM  Film  Transport 
Jodon  Engineering  Associates 

Siren 

University  Sound,  Model  1D-60T 
60  Watt,  Multi -Impedance 
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TABLE  1 

EXPERIMENTAL  - ANALYTICAL  CORRELATION  (STEEL  DISK) 


From  NASA  SP-160  "Vibration  of  Plates"  Tables  2.5,  2.29,  and  2.30 

Table  I - Experimental  - Analytical  Correlation  (Steel  Disk) 
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TABLE  2 

EXPERIMENTAL  ANALYTICAL  CORRELATION  (ALUMINUM  DISK) 


MODE 

In 

Os 

On 

Os 

2n 

Os 

3n 

Os 

4n 

Os 

On 

Is 

In 

Is 

5n 

Os 

2n 

Is 

EIGENVALUE  X FREQUENCY  HZ 


Leissa  * Experl-  Leissa  Experi- 
mental mental 


3.1 -3.5 


12.4 


21.6 


25.3 


3.3  75-85 


Dl 


I 


5.6  134-136  137 


* From  NASA  SP-160  "Vibration  of  Plates",  Tables  2.5,  2.29  and  2.30 

Table  II  - Experimental  Analytical  Correlation 
(Aluminum  Disk) 
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DISK  SPEED  n 


Figure  1.  Frequency  Speed  Diagram 
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Figure  8 - Steel  Disk 


Figure  9 - Vibration  Instrumentation 
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3N  = 3 diametral  nodal  lines 
(See  figure  1 for  explanation) 


Figure  12.  Resonant  Excitation  Frequency  as  a Function  of  Disk  Speed 
Modes  In,  2n,  and  3n 
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Figure  17.  Finite  Element  Model 


Figure  18.  Finite  Element  Resonant  Mode  Shapes  In,  2n,  and  3n 
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Figure  20.  Disk  Natural  Frequency  as  a Function  of  Speed 
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Figure  21  - Complex  Mode  Shape  of  Stationary 
Bladed  Disk 


